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SECTION  I 
INTRODUCTION 

An  item  of  primary  interest  in  the  monitoring  of  underground 
nuclear  tests  is  the  surface  wave  detectability  of  a  given  array  or  station. 

In  the  context  of  this  report  detection  levels  are  set  sufficiently  high  so  that 
for  every  detected  signal  it  is  possible  to  make  a  reasonable  estimate  of  the 
surface  wave  magnitude  from  the  Rayleigh  wave  appearing  on  a  vertical 
component  long-period  seismometer.  The  results  reported,  however,  are 
equally  applicable  to  measurements  obtained  with  a  horizontal  component 
seismometer.  The  most  accurate  method  for  estimating  detectability  is  the 
so-called  direct  method.  This  involves  processing  a  large  population  of 
events  and  plotting  the  percentage  of  detections  as  a  function  of  body  wave 
magnitude  m^  The  accumulation  of  a  statistically  meaningful  population, 
in  practice,  necessitates  data  collection  over  an  extended  period  of  time. 

An  indirect  method  for  estimating  detectability  has  also  been 
used  (Lacoss,  1969  and  Harley,  1971).  The  estimates  are  based  on  observ¬ 
ed  noise  levels  at  the  array  or  station  in  question.  This  technique  is  ad¬ 
vantageous  in  that  a  shorter  data  collection  period  is  required  to  provide 
reliable  results.  Perhaps  more  importantly,  if  one  has  some  information 
about  the  noise  levels  at  a  proposed  array  location,  this  technique  allows 
projection  of  the  performance  of  the  proposed  array.  This  permits  evalua¬ 
tion  of  the  cost  effectiveness  of  the  proposed  array. 

This  report  concerns  observed  inaccuracies  in  the  indirect 
method  as  compared  to  the  direct  method.  In  particular  it  has  been  noted 
that  direct  estimates  of  the  90%  incremental  detection  threshold  of  the 
nineteen-element  Alaskan  Long-Period  Array  (ALP A)  for  Eurasian  events 
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are  near  mb=  4.  5  for  events  in  Central  Asia  and  slightly  lower  for  Kuriles- 
Kamchatka  events  (Heiting  et  al,  1972).  Indirect  estimates  of  the  90%  level 
are  at  m^-  4.  1  for  A  =  60  and  at  mj3=  3.  7  for  A  =  35°.  It  is  believed  that 
these  errors  are  a  result  of  the  manner  in  which  estimates  made  as  a  func¬ 
tion  of  Mg  are  converted  to  estimates  as  a  function  of  mb.  In  the  past  a 
fixed  relationship  such  as  M  =  1.  59m  -3.  97  has  been  used  to  make  this  con- 

S  D 

version.  Alternatively  fixed  relationships  more  appropriate  to  the  array 

and  magnitude  range  of  interest  have  been  used.  In  either  case  this  type  of 

conversion  overlooks  the  fact  that  at  any  given  value  of  m  a  substantial 

b 

scattering  of  M  values  are  observed.  It  will  be  demonstrated  that  more 
s 

consistent  results  are  obtained  if  the  statistical  nature  of  the  M  :m  rela- 

s  b 

tionship  u  taken  into  account.  A  technique  identical  in  concept  to  that  dis¬ 
cussed  here  has  been  suggested  by  at  least  one  other  independent  investiga¬ 
tor  (Lacoss,  1971). 
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SECTION  II 


DESCRIPTION  OF  TECHNIQUE  AND  RESULTS 


The  indirect  method  as  it  has  been  applied  previously  will  be 
described  briefly.  Then  the  technique  for  incorporating  the  statistical  na¬ 
ture  of  the  relationship  will  be  presented  and  its  effects  discussed. 

A  population  of  noise  samples,  believed  to  be  representative 
of  the  ambient  noise  at  the  array  or  station,  is  accumulated.  In  the  case  of 
an  array  the  noise  samples  are  beamsteered  for  the  area  of  interest  and 
bandpass  filtered  with  a  filter  chosen  to  enhance  signals  from  the  region  of 
interest.  At  ALPA  a  beam  direction  of  340  degrees  and  an  0.  025  to  0.  055 
Hz  filter  are  used.  The  largest  noise  peak  in  each  one-half  hour  noise  sam¬ 
ple  is  then  measured.  It  is  assumed  that  if  a  signal  had  been  present  with  a 
peak-to-peak  value  at  least  three  dB  greater  than  that  of  the  largest  noise 
peak,  it  would  have  been  detected.  This  then  provides  an  ensemble  of  detect¬ 
able  signal  amplitudes  for  the  time  period  of  the  noise  samples.  This  is  con¬ 
verted  to  an  ensemble  of  detectable  surface  wave  magnitudes  using  the  formula: 


where: 


=  log  A/ T  +  1.  66  log  A 

A  is  the  peak-to-peak  value  of  the  smallest  detectable 
signal  in  millimicrons 

T  is  the  period  -  nominally  assumed  to  be  2  5  seconds 
at  ALPA 

A  is  the  epicentral  distance  in  degrees  at  which  detect¬ 
ability  is  to  be  estimated. 


The  percentage  of  times  that  a  given  could  be  detected  is  then  plotted 


versus  as  shown  in  Figure  II- 1.  These  data  were  obtained  from  a  popu¬ 
lation  of  22  noise  samples  at  ALPA  covering  the  time  period  August,  1971 
to  March,  1972  and  assume  an  epicentral  distance  of  sixty  degrees.  The 
figure  is  interpreted  as  an  estimate  of  the  incremental  surface  wave  detec¬ 
tion  probability  expressed  as  a  function  of  M^. 

It  is  customary,  however,  to  express  this  probability  as  a 

function  of  m^.  To  do  so  in  the  past,  fixed  m^  relationships  have  been 

used  such  as  M  =  1.59  m  -3. 97  (Lacoss,  1 969).  This  relationship  does 
s  b 

not  seem  to  fit  the  ALPA  data  very  well.  In  the  magnitude  range  mb=  4.  5 

to  5.  2  a  better  fit  for  the  ALPA  data  is  provided  by  the  formula  M  =  1.0m. 

r  s  b 

-1.1.  The  90%  detection  probability  of  Figure  II- 1  is  at  M^=  3.  0.  Using 
the  ALPA  conversion,  this  corresponds  to  an  m^  of  4.  1.  Direct  estimates 
of  the  ALPA  90%  detection  probability  place  it  at  about  m^  4.  5  for  events 
in  Central  Asia  (Heiting  et  al,  1972).  This  overly  optimistic  view  provided 
by  the  indirect  method  when  compared  with  the  direct  method  has  also  been 
observed  with  Norwegian  Seismic  Array  (NORSAR)  data. 

It  is  believed  that  this  discrepancy  is  caused  at  least  in  part 

bv  the  fact  that  use  of  the  fixed  M  :  m.  conversion  fails  to  account  for  the 

s  b 

statistical  nature  of  M  values  at  any  given  m  .  Observed  values  of  M  at 

S  D  S 

a  given  m^  may  lie  anywhere  within  a  range  at  least  two  units  wide.  In 

using  the  fixed  conversion  one  is  in  effect  saying,  for  example,  that  the 

curve  of  Figure  II- 1  indicates  that  events  with  =  2.  9  can  be  detected  75 

percent  of  the  time,  that  the  fixed  conversion  M  =  1. 0  m  -1.  1  shows  that 

s  b 

all  events  with  M  =  2.  9  have  m,  =  4.  0,  and  that  consequently  events  with 
s  b 

m,  =  4.  0  can  be  detected  75%  of  the  time.  In  reality,  however,  if  we  assume 
b 

that  events  with  4.  0  are  normally  distributed  about  a  mean  value  of  M^  = 

2.  9,  then  those  events  which  happen  to  have  a  larger  Mg  will  be  more  de¬ 
tectable  while  those  that  happen  to  have  a  smaller  Mg  will  be  less  detect¬ 
able.  Since  the  detection  probability  decreases  more  rapidly  as  decreases 
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FIGURE  II-l 

ALPA  SURFACE  WAVE  DETECTION  PROBABILITY 
ESTIMATED  FROM  22  NOISE  SAMPLES 


below  2.  9  than  it  increases  as  M  increases  above  2.  9»  the  net  result  is 

s 

that  the  true  detection  probability  for  m  =  4.  0  is  less  than  0.  75.  By  the 

b  J 

same  argument  detection  probabilities  for  m.  values  corresponding  to  M 

b  s 

values  below  the  point  of  inflection  of  Figure  II- 1  will  be  biased  low  if  the 
fixed  conversion  is  used.  It  is  clear  then  that  if  the  distribution  of  M 

s 

values  is  taken  into  account  in  converting  to  m^,  the  detection  probability 

curve  will  tend  to  flatten  and  the  90%  detection  level  will  occur  at  higher 

values  of  m  . 

b 

At  a  given  value  of  m^  the  probability  that  will  lie  within 

some  small  range  is  independent  of  the  probability  that  an  event  with  that 

M  can  be  detected.  Therefore  the  joint  probability  that  the  event  will  fall 
s 

in  a  particular  range  and  that  an  event  with  that  M  can  be  detected  is 

s 

simply  the  product  of  the  two  individual  probabilities.  The  correct  pro¬ 
cedure  is  to  sum  up  these  joint  probabilities  over  the  range  of  conceivable 

M  values  corresponding  to  the  given  value  of  m,  as  described  below, 
s  b 

It  is  assumed  that  at  a  given  m,  =  m  the  M  values  are  distri- 

b  s 

buted  normally  about  a  mean  value  determined  by  the  relationship  Mg  = 

1.  0  m-1.  1  .  Then  the  probability  that  the  value  of  M  lies  within  a  given 

s 

range  about  a  particular  value  given  that  mb  has  the  value  m, 

P(M.  -  e  <M  !£ M.  +  e  m  =  m),  is  the  difference  between  the  cumulative 
l  s  l  '  b 

normal  distribution  evaluated  at  the  end  points  of  this  range.  The  probability 

that  an  M  lying  in  this  range  can  be  detected  P(M.  can  be  detected)  is  ob- 
s  1 

tained  from  the  incremental  detection  probability  curve  of  Figure  II- 1. 

Finally,  the  probability  than  an  event  with  m,  =  m  can  be  detected  is  obtain- 

b 

ed  as  follows 

P(m  can  be  detected)  T"]  P(M.  -  £  <M  SM,  +  f|  m  =  m)  • 

.  l  s  l  I  b 

l 

P(M^  can  be  detected) 


II -4 


where  the  range  of  i  is  chosen  to  include  all  reasonable  values  of  M  . 

s 

This  quantity  is  computed  using  e  =  0.  025  M  units  and  for  discrete  values 

s 

of  spaced  at  increments  of  0.  05  units.  Smooth  curves  are  plotted 
through  the  computed  points  as  seen  in  Figure  II-2.  The  variable  para¬ 
meter  in  ihis  plot  a  is  the  standard  deviation  of  the  assumed  normal  dis¬ 
tribution  of  M  . 

s 

A  population  of  106  ALPA  events  with  4.  5  $m  <5.  2  have 

b 

been  partitioned  into  magnitude  bins  0.  1  m,  units  wide.  Within  each  bin 

b 

the  variance  of  was  computed.  As  a  result  of  the  relatively  small  num¬ 
ber  of  events  within  each  bin,  the  variance  differs  considerably  from  bin 
to  bin.  To  get  a  more  stable  estimate  of  the  true  variance  in  this  range  the 
variances  within  the  bins  were  averaged  giving  =  0.  18  or  a  standard 
deviation  of  0.43.  Thus  it  appears  that  the  curve  of  Figure  II— 2  most  appro¬ 
priate  to  ALPA  is  the  o  =  i). 4  curve.  This  curve  indicates  that  at  m  =  4.  5 

b 

the  probability  of  detection  is  0.  9.  This  is  in  excellent  agreement  with  the 
results  of  the  direct  method  for  events  in  Central  Asia  (Heiting  et  al,  1972). 


Since  the  assumed  M^:m^  relationship  has  a  slope  of  one,  the 

curves  of  Figure  II-2  can  be  used  for  distances  other  than  sixty  degrees. 

One  simply  shifts  them  to  the  right  by  1. 66  log  (  A/60)  units.  This  places 

the  90%  detection  probability  for  events  35  degrees  from  ALPA  at  m  =  4.  1. 

b 

Direct  estimates  for  the  Kamchatka-Kuriles  area  place  it  at  about  m,  =  4.2. 

b 

The  events  used  for  this  estimate  are  25  to  40  degrees  distant  from  ALPA. 

It  may  be  noted  that  the  340°  beam  direction  used  for  beamsteering  these 
noise  samples  is  more  appropriate  for  Central  Asia  than  for  the  Kamchatka- 
Kuriles  area.  In  practice,  however,  little  difference  in  output  noise  levels 
is  observed  between  Central  Asia  and  Kamchatka-Kuriles  beams.  Thus  it 


is  probably  not  necessary  to  form  new  beams  for  the  Kamchatka-Kuriles 
estimate. 
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SECTION  in  , 

CONCLUSIONS  ,  ! 

l 

:•  i 

* 

1  i  i 

Indirect  estimates  of  detection  probability  are  made  by  infer¬ 
ring  incremental  surface  wave  detection  probabilities  from  observed  noise 

l 

levels  at  an  array  or  station.  In  the  past  fixed  M  tm^  relationships  have: 
been  used  to  express  the  probability  as  a  function  of  m^.  It  has  been  ob¬ 
served  that  the  use  of  this  procedure  at  ALFA  and  NORSAR  leads  to  esti¬ 
mates  of  the  90%  detection  level  which  arp  about  0.'5  m^  units  lower  than 
directly  obtained  estimates  1 

i 

t 

The  use  of  a  fixed  M  :  m,  conversion  overlooks  the  fact  that 

s  b  i  i 

at  any  given  value  of  m  ,  a  wide  range  of  M  values  are  normally  observed. 

O  S 

Thus  it  is  necessary  to  use  available  information  about  the  statistical  nature 

1 

of  the  M  :m  relationship  to  express  surface  waVe  detection  probabilities 

S  D  I 

as  a  function  of  m  .  When  this  is  done  the  indirect  estimates  are  found  to 
b  1 

be  in  very  close  agreement  with  direct  estimates.  For  events  from  Central 

i 

Asia  and  from  the  Kuriles-Kamchatka  area  as  observed  at  ALPA,  the1  direct 
and  indirect  estimates  of  the  90%  detection  level  agree,  within  0.  1  m^  units. 

The  availability  of  a  reliable  indirect  method  is  a  valuable 

:  1  i 

asset.  Noise  measurements  from  Z0  or  30  noise  samples  appear  to  pro¬ 
vide  sufficient  information  to  yield  reliable  indirect  estimates.  Reliable 
direct  estimates  seem  to  require  100  or  more  events  for  a  given  region. 
Thus  it  is  possible  to  obtain  indirect  estimates  for  a  new,  or  proposed  array 
or  station  much  more  economically  and  faster  than  in  the  case  of  direct 
estimates.  This  presupposes  that  the  distribution  of  M  values  at  a  specific 
m^  are  not  strongly  dependent  on  station  location  so  that  the  0.4  standard 
deviation  reported  here  can  be  validly  applied  at  new  locations. 

i 
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